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Abstract 

We consider the structure-dependent amplitude of the decay Bg I'^l^j 
{I = e, fi) in a model based on the effective Hamiltonian for bs contain- 
ing the Wilson coefficients Cj^Cg and Cio- The form factors characterising 
the matrix elements <^7|s7^(l =F75)&|-Bs) and {j\sa^y{l =F75)&|^s) are taken 
to have the universal form fv ^ Ia ^ fr ^ fBs^BsRs/i'^E^) suggested by 
recent work in QCD, where Rg is a parameter related to the light cone wave 
function of the Bg meson. Simple expressions are obtained for the charge 
asymmetry A(x^) and the photon energy spectrum dV /dx^{x^ = 2E^/Mbs)- 
The decay rates are calculated in terms of the decay rate of Bg 77. The 
branching ratios are estimated to be Br{Bg — > e~^e~j) = 2.0 x 10^*^ and 
Br{Bg fi'^fj.^^) = 1.2 X 10~^, somewhat higher than earlier estimates. 
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1 Introduction 



The rare decay Bg — ^ /+/^7 is of interest as a probe of the effective Hamiltonian for 
the transition bs — ^ l^l^, and as a testing ground for form factors describing the 
matrix elements (7|s7^(l =F 75)b\Bs) and (7|sicr^j,(l =F Ts)^!-^^) 105 0- The branch- 
ing ratio for Bg l^l~'-f can be sizeable in comparison to the non-radiative process 
Bg —>■ , since the chiral suppression of the latter is absent in the radiative transi- 
tion. We will be concerned mainly with the structure-dependent part of the matrix 
element, since the correction due to bremsstrahlung from the external leptons is 
small and can be removed by eliminating the end-point region Si+i- ~ ^1^- (For 
related studies of radiative B decays, we refer to the papers in Ref. [0.) 

Our objective is to calculate the decay spectrum of Bg using form fac- 

tors suggested by recent work in QCD These form factors have the virtue of 
possessing a universal behaviour l/E^ for large E^, as well as a universal normaliza- 
tion. These features can be tested in measurements of B~^ fi~^i"y and Bg 77. 
We derive simple formulae for the photon energy spectrum dT/dx^, = 2E^/mB^, 
and the charge asymmetry A{x.y), defined as the difference in the probability of 
events with E^ > E_ and < E_, E± being the energies. This asymmetry is 
large over most of the x^ domain. Predictions are obtained for the branching ratios 
Br[Bg —>■ e~^e~'y) and Br{Bg /i+/i~7) which are somewhat higher than those 
estimated in previous literature 0] . 



2 Matrix Element and Differential Decay Rate 

The effective Hamiltonian for the interaction bs has the standard form 

= '^V,,V:g\cf{s^,PLb)h,l 

+ Cio{s-f^PLb)hf,^5l (1) 

Cj_ - 1 

- 2—sia^uq''{mbPR + msPL)bl-f^l j 

where Pl,r = (1 =F 75)/2 and q is the sum of the and /~ momenta. For the 
purpose of this paper, we will neglect the small g^— dependent terms in Cf^, arising 
from one-loop contributions of four-quark operators, as well as long-distance effects 
associated with cc resonances. The Wilson coefficients in Eq.(0) will be taken to 
have the constant values 

C7 = -0.315 , Cg = 4.334 , Cio = -4.624 (2) 

To obtain the amplitude for Bg ^^^~7, one requires the matrix elements 
(7|s7^(l =F75)&|-Bs) and (7|sicr^j,(l =F75)&|-Bs)- We parametrise these in the same 
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way as in Ref. 



{^{k)\rs^,b\B,{k + q)) 



e e 



e*''qPk^fv{q^)/M^ 



{l{k)\s-i^-i^h\Bs{k + q)) = -ie 



<k-q 



qk^ 



[^{k)\-sia^,q^b\Bs{k + q)) = -ee^^p^e^"" q^k'^ h^q 



'j{k)\sia^^'y5q''b\Bs{k + q)) = -ie 



e^k-q-e* ■ qk^ 



fAiq')/M, 



(3) 



The form factors fv, Ia, fx and f!p are dimensionless, and related to those of Ahev 
et al i by fy = q/Mb^Ja = I/Mb^Jt = -gi/MlJ^ = -h/Ml^. The matrix 
element for Bs — > can then be written as (neglecting terms of order rris/mh) 



-eV^V,.— 



2v^7r 



Bs 



(4) 



+ i(^l{k ■ q) - (e* ■ q)k^ (^BJ^I + B^l^l^l^ 



where 



Ml 



h 



A, = C,fv + 2Cr — 
A2 = Cio fv 

Ml , 

-Bi — Cg + 2 C^ — frp 



(5) 



q^ 



B2 = Cio /a 



(In the coefficient of C7, we have approximated ruhMB^ by Mb J. The Dalitz plot 
density in the energy variables E± is 



dT 



dE+dE^ 256tt^Mb 



(6) 



spin 
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where | 



spin 



aGi 



d^ir + l^il') 

+ (iA2r+i52p 



■ kf + (p_ ■ kf} + 2m;^(g ■ kf 
qmp+-ky + {p_^-kY}-2mKq-kf 



(7) 



+ 2Re(S*A2 + AlB2)q' 



{P+ ■ kf - ■ kf 



It is convenient to introduce dimensionless variables 

= 2EJMb, , a;± = 2E±/Mb, , A = x+ - x_ , r = m^/Mf^ (8) 

in terms of which = M^^ (1 — x^). Taking and A as the two coordinates of the 
Dahtz plot, phase space is defined by 



|A| <vx^ , v= ^1 -4mf/g2 = ^1 - 4r/(l - x^,) 
< x^ < 1 - 4r 

In terms of x^ and A, the differential decay width takes the form 



(9) 



dr 



dx^dA 



Af 



i\Al\' + \B^ 



(l-x^)(x2+A2) 1 2 

■ — rx^ 

8 2 ^ 
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(10) 



+ 2Re{BlA2 + A\B2){1 - x^)^x^A 



where J\f = [a^G^./ (25677^)] \VtbV*g\'^M^^. The last term is linear in A and produces 
an asymmetry between the and l~ energy spectra. 

We will derive from Eq. (|TOD two distributions of interest: 
(i) The charge asymmetry A{x^) defined as 



(ii) The photon energy spectrum 

+ \Bi\^){l - a; + 2r) + (l^l^ + \B2\^){1 - - 4r) 

(12) 

To proceed further, we must introduce a model for the form factors which appear 
in the functions Ai^2 and -81,2 defined in Eq. (^. 

3 Model for Form Factors 

First of all, we note that the form factors fx and /-^ defined in Eq. (|]) are necessarily 
equal, by virtue of the identity 

(y^iu = -(^fiual3(^"'^j5 (13) 

This was pointed out by Korchemsky et al . We therefore have to deal with three 
independent form factors /y, /a and /t- These have been computed in Ref. Q 
using perturbative QCD methods combined with heavy quark effective theory. For 
the vector and axial vector form factors of the radiative decay —>■ l~^i"y, and 
their tensor counterpart, defined as in Eq. (^, these authors obtain the remarkable 
result 

fyiE,) = UE,) = ME,) = ^ (q.R - ^) + (^) (14) 

where i? is a parameter related to the light-cone wave-function of the B meson, 
with an order of magnitude ~ A = Mb — rrib, where the binding energy A is 
estimated to be between 0.3 and 0.4 GeV. Applying the same reasoning to the form 
factors for Bs l^l'l, we conclude that 

ME,) = ME,) = ME,) = ^SlMjk (_g,fl, + ^)+0 (^) (15) 

In what follows, we will neglect the term Qb/^bi and approximate the form factors 

by 

f (TP \ fBs^Bs 1 1 /b 1 /-i^N 

^^'^'"^"^-^ " -iKT^s = 3 a:^; 

where A^ = M^^ — will be taken to have the nominal value 0.5 GeV. Several of 
our results will depend only on the universal form fv^A,T{E,) ~ '^/E-y, independent 



dT 



a 



7687r^ 
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of the normalization. As pointed out in a check of the behaviour fv,A ~ l/-^7 
in the case of —>■ fi'^wy is afforded by the photon energy spectrum, which is 
predicted to be 



dT 



;i7) 



In the case of the reaction Bs l^^'l, the normahzation of the tensor form factor 
friEy) at = Mb/2 {i.e. = 1) can be checked by appeal to the decay rate 
of Bs '~f'~f. To see this connection, we note that the matrix element of Bg —>■ 
•y^k, e) +7(/c', e') can be obtained from that of B^ by putting Cg = Ciq = 0, 

and replacing the factor {efTCj/q^) (J'j^l) by frix^ = l)e*'. This yields the matrix 
element 



M{Bs^^{e,kh{e',k')) 



with 

A+= -A- = ^MBjT{x, = l)Cr. 

(18) 

The result for A"^ coincides with that obtained in Refs. 0, H, H when /r(a;^ = 
1) = = lir- ^^f- B i' °f parameter is played by the 

constituent quark mass m^. ) Thus the decay width of Bg — ^ 77, 

{\A^\' + \A.\') (19) 

serves as a test of the normalization factor frix-y = 1). 

We remark, parenthetically, that the calculation of Bg —>■ '~f'~f, based on an ef- 
fective interaction for b —>■ s'-ff, produces the amplitudes A'^ and A" given in Eq. 
([T8|) in the limit of retaining only the 'reducible' diagrams related to the transition 
b — *• 57. Inclusion of 'irreducible' contributions like 6s — ^ cc — ^ 77 introduces a 
correction term in A_ causing the ratio |yl+/74_| to deviate from unity. Estimates 
in Ref. [0, H yield values for this ratio between 0.75 and 0.9. The branching ratio 
Br{Bs 77) is estimated at 5 x 10"'', with an uncertainty of about 50%. 

Having specified our model for the form factors /y(x^), /^(a;^) and /t(x^), we 
proceed to present results for the spectrum and branching ratio of Bg l'^l~'y. 
We use Mb^ = 5.3 GeV, /^^ = 200 MeV and, where necessary, A^ = 0.5 GeV in the 
normalization of the form factors in Eq. (|TBp. 
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4 Results 



4.1 Charge Asymmetry 

With the assumption of universal form factors fv = Ja = fx ^ the asymmetry 
A{x^) in Eq. (|TI]) assumes the simple form 



4 (C, + 2C,^ni -x, + 2r) + CUi - % - 4r) 

This is plotted in Fig.|l|, and is clearly large and negative over most of the domain, 
changing sign at x-y = 1 + (A negative asymmetry corresponds to /~ being more 
energetic, on average, than in the decay Bs{= bs) — > l~^l~'-f.) The average charge 
asymmetry is 

_ 3 dx,v%il-x,)2C,oiC, + 2Crj^) 

\ ) — ~A r T \ ) 

4 J^-^' dx^ vx^ {l-x^ + 2r){Cg + 2C7j^y + {1 - x^ - Ar)Cfo 

and has the numerical value (A)^ = —0.28 , (A)^ = —0.47 for the modes I = e,fi, 
the difference arising essentially from the end-point region x^ ^ 1 — Ar. 

4.2 Photon Energy Spectrum 

With the form factors of Eq. (|l^), the photon energy spectrum simplifies to 

^ = ^^^vx^!^il -x^ + 2r){C, + 2C^^^f + {l-x,- Ar)Cl^ (22) 

where the constant factor M is defined after Eq. (0). It is expedient to write this 
distribution in terms of the decay rate of 77. We then obtain the prediction 

dViB.-.in-,)ldx, \2. x3 (23) 



r(5,^77) i3 7r(l-a:^)2 



1 



{779(1 -x^) + 1} + {r/io(l -x^)} 



The first factor (in curly brackets { }) is the QED result expected if the decay 
Bs l^l^'-f is interpreted as a Dalitz pair reaction Bg 77* 1^'^^', without 
form factors. The factor (1/x^)^ results from the universal behaviour fv,A,T ~ 
given in Eq. (|TD|), while the last factor is the electroweak effect associated with the 
coefficients 779 = Cg/{2C'r) and rjiQ = Cio/(2C7). This distribution is plotted in Figs. 
and 1^, where the QED result is shown for comparison. 
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4.3 Rates and Branching Ratios 

From the photon spectrum given in Eq. (p3D, we derive the 'conversion ratios' 



^' r(B. ^ 0^) 

The numerical values are Re = 4.0% and -R^ = 2.3%. These are to be contrasted 
with the QED result given by 



MQED) = ^ 



[1 - 18r^ + 8r^) In ^ + y/T^ (-- + 13r + 4r^] 

1 — V 1 — 4r 2 



(25) 



which yields Re{QED) = 2.3%, R^{QED) = 0.67%. The absolute branching 
ratios of Bs — > ^^/~7, obtained by taking Br{Bs 77) = 5 x 10"'' 0, || are 
5r(i?, e+e-7) = 2.0 X 10-^ 5r(i?, yU+/i~7) = 1.2 x 10"^ Our resuhs for the 
average charge asymmetry (A);, the conversion ratios Ri and the branching ratios 
are summarized in Table |l|. 



5 Comments 



(i) The branching ratios calculated by us are somewhat higher than those obtained 
in previous work [|^, 0], which used a different parametrization of the form factors 
fv, Ia, fr, fx' based on QCD sum rules |jl| and light-front models [0]. In particular 
these parametrizations do not satisfy the relation Jt = fx which, as noted in [Q, 
follows from the identity cr^^ = ^e^ua/sc^^Jb- 

(ii) Our predictions for the charge asymmetry (A) and the conversion ratio T{Bs 
l^l^'-f)/r{Bs '~f'~f) are independent of the parameter which appears in the form 
factor in Eq. (|16|). The branching ratios in Table |^ assume Br{Bs '~f'~f) = 5 x 10"^, 
and can be rescaled when data on this channel are available. 



(iii) A full analysis of the decay Bg l^l^'j requires inclusion of the bremsstrahlung 
amplitude corresponding to photon emission from the leptons in Bg — > . This 
contribution is proportional to Jbs^i ^i^nd affects the photon energy spectrum in 
the small region. We have calculated the corrected spectrum for Bg l^l^l, 
following the procedure in |TT[, and the result is shown in Fig. H for the case / = /i. 



As anticipated, the correction is limited to small x^, and can be removed by a cut 
at small photon energies. 
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(iv) The QCD form factors in Eq. (p!6| ) are valid up to corrections of order (Aqcd / E.y)'^. 
In the small region, arguments based on heavy hadron chiral perturbation the- 
ory suggest form factors dominated by the B* pole with the appropriate quantum 
numbers, for example, 

^^■<^^) ~ Ml(l-l,)-Ml^ P«> 

Defining Mb* — Mb^ = AM, this form factor has the behaviour /y(x^) ~ '^;r+S' with 
5 ~ 2AM/Mb^ ~ 0.02. We have investigated the effect of replacing the QCD form 
factor of Eq. (0) by a different universal form fv,A,T{xy) = /bJ {SAs{x^ + 6)), and 
found only minor changes in the numbers given in Table |I|. In general, one must 
expect some distortion in the spectrum at low Xj, compared to that shown in Figs. 



(v) We will examine separately the predictions for A{x^) and dT/dx^ in the reaction 
Bs — > r"'"r~7, in which the bremsstrahlung part of the matrix element plays a 



significant role We will consider also refinements due to the g — dependent 



term in Cf^, and the effects of cc resonances. 

In view of their clear signature, non-negligible branching ratios and interesting 
dynamics, the decays Bg l^l^^ could form an attractive domain of study at future 
hadron colliders. 
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Decay 


Average Charge 
Asymmetry 

(A) 


Conversion Ratio 

r{B,^i+i-^) 

r(B,-^77) 


Branching Ratio 

r{Bs^i+i-^) 

r(Bs-^all) 


Bs — >■ e+e~7 


-0.28 


4.0% 


2.0 X 10-s 


Bs //+/i~7 


-0.47 


2.3% 


1.2 X 10-^ 



Table 1: Average charge asymmetry, Conversion ratio and Branching ratio for the 
decays Bs — > e'^e"^ and Bg — > (Last column assumes Br{Bs 77) = 

5 X 10-7) 
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Figure 2: Photon Energy Distribution for Bg e"'"e~7, normalized to Bg 
(Dashed hue is the QED result.) 
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Figure 3: Photon Energy Distribution for Bg /i+/i^7, normalized to B, 
(Dashed line is the QED result). 
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Figure 4: Photon Energy Spectrum in Bg — > /I'^iJ, 7, with bremsstrahlung (sohd 
hne) and without bremsstrahlung (dashed hne) 
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